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Abstract 

Several modeling approaches can benefit from data and dynamics transformations that allow simplification, general¬ 
ization, and/or more conceptual treatment of the studied problem. In this work, we discuss the phasor representation, 
capable of mapping sinusoidal signals into values in the complex plane, which can be associated to vectors. As such, 
phasors allow visualizating and understanding amplitude and phase relationships, such as relative phase and phase 
coherence. The potential of phasors for diverse applications is illustrated with respect to AC linear circuits, filters, 
oscillators and the Kuramoto model. 


“...this brave o’erhanging firmament, this majestical roof 
fretted with golden fire.’ 

William Shakespeare. 


1 Introduction 

To a good extent, scientific modeling as well as technolog¬ 
ical developments often benefit from approaches capable 
of simplifying, and/or generalizing and/or obtaining more 
effective and approachable representations of aspects of 
the studied phenomena. This is frequently the case with 
transformations such as the Fourier series, in which a wide 
range of signals can be expressed in standardize way, more 
specifically in terms of a sum of sines and cosines with spe¬ 
cific amplitudes, frequencies, and phases. Other examples 
include the Laplace transform, the Taylor series, the use 
logarithm to map products into additions, and the use of 
moments to express statistical properties of signals. 

In the present text, we focus on phasor representations 
of signals. Basically, this approach allows a sinusoidal 
signal of time - with amplitude A , frequency / (or the 
associated angular velocity uj = 2nf) and phase 0 - to 
be mapped into a respective complex number, the pha¬ 
sor <£> = Ae l °, where i = T (e.g. [I]). Recall that this 
complex number can also be associated to a vector in the 
plane. Phasor representations disregard explicit represen¬ 
tation of time as well as the signal angular velocity. These 
two quantities can nevertheless be taken into account by 


assuming that the mapped phasors are rotating around 
the complex plane origin with angular velocity uj. The 
concept of phasor is illustrated in Figure [l] 

The advantage of the phasor approach derives mostly 
from its ability to represent in an intuitive and simple way 
the phases of the signals, as well as their interrelationships 
(e.g. relative phase, phase coherence, etc.). In addition, in 
cases where all signals have fixed frequency, it is possible 
to perform several operations between the involved signals 
in an objective and simplified geometrical way (e.g. vector 
addition). 

As a consequence of these interesting properties, pha¬ 
sors are adopted and used in a wide range of scientific 
and technological areas, including but by no means be¬ 
ing limited to: physics (e.g. oscillations), dynamical sys¬ 
tem, electrical and electronic engineering, dynamical sys¬ 
tems, control theory, complex systems, vibrations, acous¬ 
tics, audio, neurosciences, among many others. 




Figure 1: Plot of the signal x(t) = l ot -\- 6 (a), and its respective 
phasor representation in the complex plane (b). 

We shall start by revising some important properties 
of the sine and cosine functions, introducing the concept 
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of phasor representation, and then illustrating its poten¬ 
tial with respect to some specific applications, namely AC 
(alternating current) linear circuits, linear filters, individ¬ 
ual harmonic oscillators, systems of harmonic oscillators, 
and the Kuramoto model (e.g. mm- The latter corre¬ 
sponds to a particularly interesting approach to modeling 
and understanding some remarkable synchronization phe¬ 
nomena found in nature including fireflies flashing, heart 
and neuronal cell operation, clocks synchronization, and 
even aspects of economical and social dynamics. 

2 The Sine and Cosine Functions 

The sine and cosine trigonometric functions are particu¬ 
larly important in science, as they model several aspects 
of the real world, such as several types of oscillations. 
They are also frequently used in several technological ap¬ 
plications. Though intrinsically simple, these functions 
exhibit some interesting, more subtle, features often aris¬ 
ing from their periodical behavior. Some of these prop¬ 
erties, namely those more closely related to phasors, are 
discussed in this section. 

Both the sine and cosine functions take an angle (mea¬ 
sured in radians or degrees) as argument, which we will 
call 4>(£), so that we have sin(<$>) and cos(<h), respec¬ 
tively. Their 27r-periodicity can now be expressed as 
sm(<h) = «sm(<h + 27r) and sin($) = sin(^-\-2 tt). Though 
<F can take any real value, it is often interesting to con¬ 
sider that 0 < < 27r. 

One of the first interesting properties to be noticed 
about the sine and cosine functions is that they are closely 
interrelated, in the sense that any of them can be obtained 
by displacing the other along the angle argument. For in¬ 
stance, we have that 


sin(&) = cos(& — tt/2), (1) 

cos(<F) = sin(& + 7t/2). 

It follows from this property that what we discuss and 
learn about one of these functions can be immediate and 
easily transferred for the other. In this work, we will focus 
on the sine function without loss of generality. 

Several of the applications of the sine function involve 
expressing the angle argument as a function of time, 
i.e. <F(£). Typically, this time-dependent argument is ex¬ 
pressed in terms of a respective frequency f and a phase 
0. The period of this function is T = 1//. Observe that 
the phase 0 corresponds to the value of the sine argument 
<!>(£) at the initial time, i.e. <f>(£ = 0) = 6. The angular 
velocity (or frequency ) uj = 2irf is also often considered. 
So, we can define a sinusoidal oscillating signal xft ) as the 
following function of t: 


xft ) = A sin (2 tt ft + 0) = A sin (cot + 6) . (2) 

One example of this type of signal xft) is illustrated in 
Figure [lja). 

While the frequency / determines how fast the sig¬ 
nal changes, the phase 0 can be understood as a con¬ 
sequence of displacements with respect to the reference 
sinusoidal signal sin(27rf). The frequency / (or the an¬ 
gular frequency uj) and the phase 0 can be understood 
as representing parameters of the function in Equation [2] 
Unlike variables , such as the independent-value t or the 
dependent-value xft), parameters remain constant during 
a given realization of the function. For instance, a lin¬ 
ear pendulum will oscillate while keeping its parameters 
/ and 6 constant. 

We have that / is commonly measured in Hertz ( Hz ), 
uj in rad/s (radians per second), and 0 in rad. Note that 
degrees (°) can also be used instead of radians. 

An interesting situation is defined when the phase is a 
function of time, i.e. 

xft) = A sinfjjt + Oft)) (3) 

In particular, if Oft) is an affine function of t, i.e. Oft) = 
vt + 0q , we have 


xft) = A sinfjjt + vt + 0q) = A sin(uj e t + 0q) (4) 


with uj e = [uj + v\. In other words, in this situation 
the phase Oft) behaves in an analogous way to an angular 
frequency. As an example, consider the situation in which 
we want to obtain x^ ft) = sin(2ujt) from a signal x\ = 
sinfjjt + Oft)). This can be easily achieved by making 

v — UJ. 

Recall that we have referred to the argument of the sine 
function as <£>(£). Now, the respective time derivative of 
this argument, given as 


d$ft) 

dt 


= m 


(5) 


has unit rad/s and can be effectively understood 
as the instantaneous (angular) frequency of the signal 
Asin(&(t)). For instance, if Oft) = vt + we have, as 
expected, that 4>(t) = d(ujt + vt + 0q)/ dt = uj + v = w e . 


3 Phasors 

It is often interesting and useful to visualize signals such 
xft) in Equation [l8|a). We have already illustrated 
the Cartesian representation of the sine function in Fig¬ 
ure [lja). However, this type of representation can be 
highly redundant as a consequence of the 2tt periodicity 
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of the sine. A more effective approach can be obtained by 
using the Euler equation to map the sine function in the 
complex plane (e.g. imimiT]), i.e. 

x(t) = Asin($>(t)) = Aim (6) 

which can be associated to the time-varying complex 
variable z(t) = Ae % ^\ 

Figure [ljb) shows the representation of the signal 
x(t) = Asin{2^ft + 0) in Figure [lja) into the complex 
plane. Observe that the sine values are mapped into a 
circle of radius A , while its argument 0 < <L(£) < 2 tt cor¬ 
responds to the angle between each of the function points 
and the Real axis. It is also possible to understand the 
respective signal x(t) as being obtained by rotating the 
phasors counterclockwise with constant angular velocity 
uj while taking the respective projections of the phasor 
onto the imaginary axis. 

For the very same reason that the phasor representation 
allows the elimination of the redundancy derived from 
the sine periodicity, it also implies that we are no longer 
able to geometrically discriminate between values of <L(£) 
(i.e. angles) that are larger than 2 i r or smalller than 0. 
In other words, we can only visualize the function x(t) 
as defined for 0 < <£>(£) < 2 i r. The discussed Cartesian 
and complex representations of x(t) can be understood 
as being unwrapped and wrapped around the period 0 < 
$>(t) < 2 i r, respectively. 

Now, the smallest difference between the two phases, 
which provides a measurement of the similarity between 
the two respective phases, can be immediately appreci¬ 
ated in terms of the smallest angle between the two pha¬ 
sors in the respective complex-plane representation. Note 
that, in the general case, this angle cannot be always ob¬ 
tained mathematically by making S = #2 — 0 \, as a conse¬ 
quence of the periodicity of the phasor representation. For 
instance, in case 0\ — 5° and 62 = 355°, we would have 
8 = 350° instead of the expected smallesr angle differ¬ 
ence of just 10°. For such reasons, it becomes interesting 
to use other ways of obtaining measurements expressing 
angle differences while disregarding the possible discon¬ 
tinuities implied by the angle periodicity. An interesting 
possibility is to adopt A 0 = \sin {02 — 0 1 )|. 

4 AC Linear Circuits 

A powerful application of phasors regards the time- 
varying analysis of linear electric circuits, which is usually 
called alternated current (AC) analysis (e.g. [11 8]), as a 
counterpoint to direct current (DC) analysis. The AC 
regime involves a transient period followed by an steady- 
state stage. Here we focus attention on the latter situa¬ 
tion, characterized by the almost complete elimination of 


Table 1: The reactances and impedances of the passive components 
R , C and L. Observe that i = ed 7 */ 2 ) and —i = ed -7r / 2 ). 


component 

reactance 

impedance 

R 

R 

R 

C 

X C = ~^c 

z c = iX c = ~zjc 

L 

X L = luL 

Zl — iXL = iuL 


the effects of the initial charges and fields in the capaci¬ 
tors and inductors. The complete behavior of the circuit 
can be obtained by linear superimposing the transient and 
steady-state solutions. 

Interestingly, the linear circuit framework is not only 
very useful in electricity and electronics, but it can also 
be related and applied to other concepts and areas such as 
acoustics, vibrations, temperature, flow, hydrodynamics, 
magnetism, among others. 

AC linear circuits involve combinations of some of the 
following possible components: voltage sources (v(t), in 
Volts), current sources (i(£), in amperes, A), resistors 
(R, in ohms, Q), capacitors (C, in farads, F), and induc¬ 
tors (L, in henries, H). Sources are understood as active 
components, while resistors, capacitors and inductors are 
said to be passive. Note that, in addition to at least one 
time-varying voltage or current source, AC analysis may 
involve also other time-varying or time-constant voltage 
and/or current sources. Circuits incorporating more than 
one source can be analyzed through linear superimposition 
of the solutions found considering each of them separately. 

The time-varying signals are assumed to be in the har¬ 
monic form. For instance, we can have a voltage source 
v(t) = A sin((dt + 9), where A G 5ft is the signal amplitude. 
This signal can be represented in the complex plane as 

v(t) = Im (a e i(a;t+0) ) = Aim ( e iujt . e id ) (7) 

So, the original voltage signal v(t) is transformed into 
the real number A. e luJt . e 10 . Because a linear circuit with 
a source operating at uj will have currents and voltages at 
all passive components oscillating with that same angular 
frequency, we can take into account only the amplitude 
and phase of these signals, transforming v(t) into the re¬ 
spective phasor V = Ac 10 . 

In a purely resistive resistor R, the voltage and current 
across it can be related by Ohm’s law v(t) = Ri(t). This 
law can be immediately extended to capacitor and induc¬ 
tors by defining the respective reactances and impedances , 
as given in Table [l] As with resistance, the reactances and 
impedances are measure in Ohms. 

Observe that reactances assume real values, while 
impedances are purely imaginary. Let’s consider a com¬ 
ponent consisting of a resistor R in series with a capacitor 
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C and an inductor L. The total reactance is defined as: 


(a) 


X = Xc + X L (8) 

and the total impedance is given as 

Z = R + iX = R + i(Xc + X L ) = i2 + T Zfi. (9) 

This composed impedance device represents a general 
case of impedance in a linear AC circuit. Interestingly, 
impedances Z can be associated in a similar way to finding 
the equivalent resistance of a network of resistors. 

The above definitions allows the Ohm’s law to be 
rewritten as: 


V = ZI (10) 

in which / is oriented in opposition to V along the pas¬ 
sive component Z. 

Let’s consider a single capacitor C connected to a volt¬ 
age source V = Ae l ° with a generic phase 0 and ar¬ 
bitrary amplitude A. For this case, we have that I = 
V/Z = VjZc = V/(-i/uC ) = iVuoC = e^^e^AuC = 
cjjACe h 7r / 2 + 6 9, an d consequently the voltage V lags (in the 
sense of being delayed in time) the current / by a relative 
phase 7r/2, irrespectively of cj, A , and C. The respective 
phasors are illustrated in Figure [2|a). 

Now, if instead of a capacitor, we attach an inductor 
L to the source V, we have I = V/Zl = V/(iujL) = 
—iV/(<jjL) = Ae^-^e 0 /(ljL) = A/(uL) e^""/ 2 ), and 
now it is the current that lags the voltage by a relative 
phase 7r/2. The respective phasor diagram is given in 
Figur^b). 

As a more general example, consider the circuit con¬ 
taining a resistor R = 100011 and a capacitor C = 0.01 ptF 
in series with a voltage source v(t) = 2sin(2irft + 7r/4), 
i.e. V = 2e*( 7r / 4 ) Volts , as shown in Figure^ We adopt 
/ = 1/(27 tRC) ~ 15915.49RT, so that c<j = 10 5 rad/s, 

It follows that Xc = —l/(10 _3 )fl = —10 3 fl and Zc = 
-i 10 3 fl, so that Z = lOOO-ilO 3 ^ = 1000+10 3 e^- 7r / 2 )fl. 
Now, we can determine the current I by using the complex 
Ohm’s law as: 


I = 


V 

Z 


2 e hV 4 ) 

1000 + 10 3 ed-^/2) 


V2e i(7r/2) = i V2 mA. 


( 11 ) 

By considering R and C in terms of their separated 
impedances Z R = 7? = 100011 and Zc = —i/(uC) = 
-il0 3 S2, we can now obtain the voltages across these com¬ 
ponents as Vr = RI = y/2e 1 ^^ volts and Vc = Zcl = 
10 3 a/2 10 -3 = y/2 volts. Observe that we could also have 
calculated Vc from V and Vr by using Kirchhoff’s voltage 
law, i.e. V = Vc + V R ^V c = V-V r = 2e < ( 7r / 4 ) -i V2 = 
\/2 volts. 



Figure 2: (a): In a capacitor , the phasor notation helps to visualize 
that the voltage V lags the current / by 7r/2. (b): In an inductor , 
the current / lags the voltage V by n/2. (c): The voltage phasors 
with respect to theexample in the text, consisting of voltage source 
V in series with a resistance R (voltage Vc) and a capacitor C 
(voltage Vr). The current phasor /, not shown, is parallel to Vr. 


Note that Kirchhoff’s voltage law can be suitably rep¬ 
resented as the complex (or vector) sum in Figur^2|c). 
Observe also that the phasor I defined by the current i{t) 
is parallel to the phasor associated to Vr. 


5 Linear Filters 

AC circuits such as that in the RC example can be under¬ 
stood as filters (e.g. ID 0), in the sense that they will at¬ 
tenuate the input signal (voltage or current) as a function 
of the signal frequencies uj. Filters are very important 
in several applications ranging from electronics to neu¬ 
rosciences. In this section we illustrate how the phasor 
approach can be applied to better represent and under¬ 
stand filter operation. Generally speaking, linear filters 
involve the passive components resistors, capacitors and 
inductors. A filter is often studied by considering its op¬ 
eration, i.e. how it changes the magnitude and phase of 
the input signal V in terms of the respective frequency /, 
which can be expressed as: 


where ay is the voltage gain and cj)y is the relative 
phase implied by the filter. Observe that both these gains 
are functions of the frequency /. 

An example of a linear filter corresponds to the series 
configuration of V, R and C addressed in section @ By 
allowing V to vary in terms of /, instead of fixating it 
at / = l/(27rRC), we have that the impedance of the 
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capacitor will also change with / as Zq = —i/(u)C), so 
that / = V/Z = V/{R - i/(wC)), implying V c = Z C I = 
V/(l + iujRC ), so that 


Vc = 1 

y 1 + iujRC 


(13) 


Figure [3] shows the phasors V, Vr and Vc as the 
frequency / increases, in steps, from /i = 500 Hz to 
fg = 50 kHz along the sense indicated by the respec¬ 
tive integer indices. Observe that, for any frequency, we 
can obtain the phasor associated to the voltage source as 
V = V R + Vc, as illustrated by the dashed black arrows 
in this figure for fg. It can be readily verified from the 
obtained phasor representation that, for small values of 
f, V R tends to align with V ; while as / increases, Vc 
will be substantially reduced while aligning with V as a 
consequence of the reduction of the capacitive reactance 
verified for higher frequencies. 



Figure 3: The behavior of the RC filter example from Section^](R = 
1000Q, C = 0.01 rF and V = 2e Z7r / 4 can be intuitively represented 
in phasor notation. 


As the magnitude of Vc attenuates with /, it becomes 
interesting to define some critical point attenuation point. 
The standard reference value is the cutoff frequency, f c , 
corresponding to where the signal power is halved. This 
implies an attenuation by ay = y/2/2 ~ 0.707 of the 
signal. It can be found that 


fc = 


(14) 


2nRC 

For f = f c , we also have that the signal corresponding 
to V undergoes an attenuation of 20Lo^(\/2/2) ~ —3 dB 
and the output signal lags V by —7r/4 (as is the case with 


the previous example, see Figure |2|. The passing band 
of this RC filter can the be understood as the frequency 
interval 0 < / < f c , implying this circuit to operate as a 
lowpass filter. 


6 Harmonic Oscillators 

Let’s now briefly consider an important physical phe¬ 
nomenon, namely harmonic oscillations. Consider the 
pendulum (e.g. 1) in Figure[4ja). Here, a sphere (or bob) 
with mass m is attached to a pivot through a weightless 
rod of length L which can rotate freely around its pivot. 
For simplicity’s sake, we will assume that all the mass m 
is concentrated at the center of mass of the sphere. Also, 
it is assumed that the pendulum is in the vacuum and 
that the pivot is frictionless, so that there is no energy 
dissipation. The angular position /3(t) of the sphere is 
taken with respect to the vertical, and the initial angular 
position is taken as /3(0) ■= /3g. The initial angular veloc¬ 
ity is assumed as /J(0) = 0, i.e. the sphere is initially at 
rest. 



Figure 4: Modeling the angular movement of a pendulum. A sphere 
with mass m is attached to a pivot through a weightless rod. The 
angular position (3{t) of the sphere is taken with respect to the 
horizontal axis. 


The mass is under the force of gravity / = mg, which 
can be decomposed into two orthogonal components, one 
along the rod orientation, f n , and another tangential to 
the sphere displacement, f t . The action of f n is counter¬ 
acted by a respective reaction f r , as there is no possibility 
of the rod to expand or contract, so that f n + f r = 0 and 
then f r = —f n . On the other hand, the tangential com¬ 
ponent induces an angular displacement of the sphere. 
By trigonometric decomposition of /, we derive the mag¬ 
nitudes f t = f sin(f3(t)) and f n = f cos(/3(t)). Recall 
that the angular acceleration j3(t) is related to the linear 
counterpart a(t) as a(t) = L/3(t). By Applying Newton’s 
second law, we derive: 
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ft = mg sin(/3(t )) = —m a(t ) ~ —m L f3(t) => (15) 

Pit) = -sin{P{t)) 

Though this differential equation completely specifies 
the angular displacement /3(t), it is difficult to solve as 
a consequence of its nonlinearity. For small angular dis¬ 
placements /3(t) <C 1 rad, we can adopt the approxima¬ 
tion sin(/3) ~ /3, implying the simplified linear differential 
model 




(16) 


We observe that the name harmonic motion tradition¬ 
ally refers to the fact that the restoring force is directly 
proportional to the opposite displacement, as in the pre¬ 
vious equation. 

yields /3(t) = 


The solution of Equation 16 


/3qcos {^^/g/Ltj = /?o sin \^Jg/L t + 7r/2^. Thus, 

we find that the considered pendulum undergoes a sinu¬ 
soidal oscillation around the vertical axis with angular 
amplitude of the oscillation equal to /3 q. By considering 
that yfgjL = c a(t), we have / = 1/(27 x)^/gjL. By 
making 0 = 7r/2 , the solution can now be expressed as 


/3(t) = /3q sin(ut + 6) (17) 

which is in the same form as Equation [2] 

Therefore, we obtain that a linear pendulum produces 
a signal that can be described by a phasor. This is partic¬ 
ularly convenient in the case of systems involving several 
oscillators, as we will see in the next sections. 


7 Independent Harmonic Oscilla¬ 
tors 

Let’s now consider not one, but N oscillators that are 
completely disconnected, and therefore, independent one 
another. Each of these oscillators j will have respective 
angular frequency ujj and phase Oj. Because the phase 
is 27r-periodical, it is convenient to visualize the phase in 
the complex plane. Considering the Euler equation, we 
have: 

sin (ujt + Oj) = Im (e i ^ t+e ^ = Im (e i( ^ 4) . e^) 

(18) 

In case all the angular frequencies are equal , meaning 
that ujj = uj for any oscillator j, the phase difference be¬ 
tween any pair of oscillating signals Xi(t) is time-invariant, 
so we can conveniently visualize the oscillations only in 
terms of the respective quantities e^'(i.e. a phasor) de¬ 
fined by the respective phases Oj. 


Figure [5] shows, in the complex plane, the distribution 
of N — 100 phases following a uniform distribution in the 
interval [0,27r). In order not to clutter the visualization 
of the phasors, each of them is represented only in terms 
of small circles of different collors. 



Figure 5: Phasor representation of the phases of N = 100 oscil¬ 
lators of the type Xj{t ) = sin(ujt + 0j), where the phases 0i were 
drawn with uniformly random probability in the interval [0, 27t). 
The black dot near the center of the complex plane corresponds to 
the parameter II in Equation |19| Note that n —> 0 as N —* 00 . 


This type of phase distribution presents an intrinsic ‘or¬ 
der’ that can be evaluated in terms of the parameter: 


n = re^ = 


1 

N 


N 

3 =1 


(19) 


where r gives the phase coherence and is the average 
phase. Observe also that II corresponds to the phasors 
center of mass. For uniformly random distributions such 
as that in Figure [5j it can be shown that r —>> 0 as N —>> 00 . 


8 The Kuramoto Model 

The Kuramoto-Daido model, or simply Kuramoto model , 
was developed by Yoshiki Kuramoto along the 70’s as 
a simplified version of Winfree’s synchronization ap¬ 
proach [2j. Kuramoto was then able to solve this non¬ 
linear dynamical problem in analytical fashion. Let’s con¬ 
sider a set of N coupled oscillators exhibiting dynamics 
in the form 


Xi(t) = sin (uit + 0i(t)) . (20) 

Observe that the phases can vary along time. We 
will express the argument of each oscillating signal as 
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®i(t) = oJi + Oi(t). In the Kuramoto model (e.g. [3] 0]), 
the dynamics performed by each oscillators is assumed to 
be governed by the following equation 


N 

9i(t) = ui + K id sin (6j{t) - O^t)) (21) 

2=1 

for i = 1, ..., TV. Here, if is a coupling matrix and 
uji and Oi(t) are, respectively, the angular velocity and 
phase of each oscillator i. Observe that when K = 0, 
i.e. the oscillators are uncoupled, we have 6 i = c^, and 
each oscillator i will acquire an effective angular velocity 
2 ui. 

Let’s take into account the situation arising when each 
oscillator has a slightly different frequency More 

specifically, it is assumed that these values are drawn 
from a normal distribution with average (i and relatively 
small standard deviation a. For simplicity’s sake, but 
without loss of generality, we make // = 0. Let’s also 
assume that all coupling weights are equal to a value 
if /TV, i.e. Kij = K/N , and that if is relatively small, 
e.g. Km 1. As all oscillators in the Kuramoto model have 
the same amplitude, all phasors in the complex plane will 
belong to the circle of radius 1 centered at the origin of 
the plane. 

Figure [6] illustrates the unfolding, along time, of the 
phase distribution (calculated numerically) represented in 
terms of the respective phasors. Here, K = 1, /i = 0 and 
a = 0.1 rad. 

Recall that the angular velocities Ui of each oscillator 
are slightly different one another, remaining fixed along 
time. So, how can these osccilators synchronize in fre¬ 
quency? This takes place because of the phase variations 
implied by the governing Equation [2l] effectively act on 
the respective argument <F(£), i.e. uit + 6 i(t), compen¬ 
sating the respective angular velocity 0 i so as to achieve 
similar resulting angular velocity aj e ^(t) = d$(t)/dt = 
uji + dOi (t) /dt (see Section [2| . 

As a consequence, even if the oscillator phases are not 
strongly clustered, the effective oscillator frequencies can 
be synchronous. When the phases tend to merge along 
time, as in Figure [6j we have the phenomenon of phase 
locking , which is favored by higher coupling values. 

An important aspect of the Kuramoto dynamics re¬ 
gards its specific stability. Depending on r and if, the syn¬ 
chronization can be stable or unstable , defining respective 
bifurcation diagrams (e.g. mm- The Kuramoto model 
has motivated several extensions and variations, including 
the consideration of complex network topologies (e.g. (TO]) 
on the respective dynamics (e-g. 0). 


9 Concluding Remarks 

Despite their relative simplicity, phasors allows us to vi¬ 
sualize and identify relationships between phases in sev¬ 
eral types of problems, including circuits, oscillators and 
complex systems. In this text, we presented the phasor 
notation at an introductory level and illustrated its po¬ 
tential for applications with regard to an RC circuit, the 
harmonic oscillations of a linear pendulum, a system of in¬ 
dependent oscillators, and the Kuramoto model. Several 
other applications exist, and others will be devised. 
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Costa’s Didactic Texts - CDTs 


This is a Costa’s Didactic Text (CDT). CDTs 
intend to be a halfway point between a formal sci¬ 
entific article and a dissemination text in the sense 
that they: (i) explain and illustrate concepts in a 
more informal, graphical and accessible way than 
the typical scientific article; and, at the same time, 
(ii) provide more in-depth mathematical develop¬ 
ments than a more traditional dissemination work. 

It is hoped that CDTs can also provide integration 
and new insights and analogies concerning the 
reported concepts and methods. We hope these 
characteristics will contribute to making CDTs 
interesting both to beginners as well as to more 
senior researchers. 

Though CDTs are intended primarily for those 
who have some preliminary experience in the 
covered concepts, they can also be useful as 
summary of main topics and concepts to be learnt 
by other readers interested in the respective CDT 
theme. 

Each CDT focuses on a few interrelated concepts. 
Though attempting to be relatively self-contained, 
CDTs also aim at being relatively short. Links 
to related material are provided in order to 
complement the covered subjects. 

The complete set of CDTs can be found 
at: https://www.researchgate.net/project/ 

Costas-Didactic-Texts-CDTs. 
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Figure 6: Time unfolding of phasors obtained from the phases of N oscillators in the Kuramoto model assuming K = 1, /i = 0 and a = 0.1 
rad. The angular velocities were drawn from a normal distribution with null mean and unit variance, and are remain fixed along time. The 
plot titles indicate the respective time interaction, while the red X identifies the phasors center of mass, leaving the trajectory in black. 
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